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Is^ ■ Abstract. The '^He('^He, 7)^Be and ^H(^He,7)^Li reactions are studied in an 

Q"^ [ extended two-cluster model which contains a + h/t and ^Li -|- p/n clusteriza- 

^N I tions. We show that the inclusion of the ^Li + p/n channels can significantly 

^S^ ■ change the zero-energy reaction cross sections, S'(O), and other properties of 

I , the ^Be and ^Li nuclei, like the quadrupole moments Q. However, the results 

CJ I agree with the known correlation trend between ^(O) and Q. Moreover, we 

H ■ demonstrate that the value of the zero-energy derivatives of the astrophysical 

• • , S-factors are more uncertain than currently believed. 

Vh ' 1 Introduction 



The ^H(^He, 7)^Li and ^HeC^He, 7)'^Be reactions play important roles in astro- 
physics. The former process is an important ingredient in big-bang nucleosyn- 
thesis, while the latter one is a key reaction in the solar p-p (proton-proton) 
chains. The ^He(^He,7)^Be reaction produces ^Be in our Sun, which is destroyed 
by either the capture of an electron or a proton in the ^Be(p, 7)^B process. In 
both cases, neutrinos are generated which are observable by terrestrial detec- 
tors and thus contribute to the solar neutrino problems [y. The production rate 
of ^Be in the Sun is determined by the competition of the '^He(^He,7)^Be and 
^He(^He, 2p)^He reactions. Therefore, the precise knowledge of these cross sec- 
tions at the most effective solar energies (~ 20 keV) is an important input to solar 
models and to the possible solution of the solar neutrino problems. 

There exist several measurements on the ^He(^He,7)^Be cross section down 
to energies of about 100 keV. From these data, the reaction rate at solar en- 
ergies is obtained by extrapolation. The ^He -|- ^He fusion reaction at low en- 
ergies is usually viewed as an approximate external capture process |2[. Thus, 
*Dedicated to Achim Weiguny on the occasion of his 65th birthday 



2 Study of the ^He(*He, 7) ''Be and ^H(''He, 7)'^Li reactions 

the energy dependence of the cross section can then be obtained from a simple 
^He + ^He potential model, which in fact reproduces the energy dependence of the 
measured data rather well. The most recent compilation of the data quotes the 
^He(^He, 7)^Be cross section as 5*34(0) = 0.53 ± 0.05 keVb |y], where the rather 
small uncertainty mainly reflects differences in the absolute normalization of the 
cross section. As customary, we have used the standard S-factor parametrization 
of the cross section (the subscript '34' stands for '^He + ^He) 

S{E) = aiE)Eexp [27rr?(i?)] , rj{E) = ^^^f . (1) 

where Zi and Z2 are the charge numbers of the fragments in the incoming channel, 
while ^ and k are their reduced mass and wave number, respectively. 

There exist also more elaborate theoretical studies of the ^H(^He, 7)^Li and 
^He(^He,7)^Be reactions in the literature, which are based on the resonating 
group method or the related microscopic potential model |y, ^, ^, Q |8|, |§, |lO|] . 
These microscopic studies, based on a single ^He + ^He fragmentation, yield an 
energy dependence of the ^He(^He,7) Be cross section, which is similar to the 
potential models. Furthermore as a success of these approaches, they correctly 
predicted |^, ^ the low-energy ^H(^He, 7)^Li cross section in absolute normal- 
ization and energy dependence, which has been subsequently confirmed by the 
precision measurement of Brune et al. |]lT| . 

Why do we then want to study these low-energy reactions again? Despite its 
apparently successful agreement with data, the dependence of the cross section on 
the enlargement of the model space has not been yet established. A first attempt 
has been made by Mertelmeier and Hofmann |^] who included additionally a 
®Li -|- p configuration in their model space. This multichannel RGM calculation 
shows a somewhat increased S-factor compared to the single-channel case, while 
the rise towards lower energies is reduced. Given the astrophysical relevance it 
is important to determine whether this weaker energy dependence is physical or 
an artifact caused by the approximate treatment of the asymptotic Whittaker 



functions in the bound states. Second, we have found in Ref. [12| that taking into 
account the ^Li -|- p clusterization changes some of the key properties of the Be 
bound states considerably. For example, the ^Be quadrupole moment is increased 
by about 0.5 — 1.0 efm^ if the ^Li -|- p configuration is added to the dominant 
^He -|- ^He clusterization. 

In the present work we study the ^He(^He, 7)^Be and '^H(^He, 7)^Li reactions 
in an extended two-cluster model, which can take into account the ^Li -|- p/n 
clusterization, in addition to the more important ^He -|- ^He/^H configurations. 
Our main goal is to explore how much the energy dependence of the low-energy 
^He(^He, 7)^Be cross section is affected by the addition of the ^Li + p configura- 
tion. The simultaneous study of the ^H(^He, 7)^Li reaction serves as a potential 
check, as in this case the low-energy cross section is known rather well. 
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2 Model 

Our study of the ^He(^He, 7)^Be and ^H(^He, 7)^Li reactions is based on the mi- 
croscopic cluster model. As a model space we adopt the dominant ^He + ^He/^H 
clusterization, to which the ^Li+p/n configuration is added as well. For compari- 
son, we also perform calculations based solely on the ^He-|-^He/^H clusterization; 
in this case our studies correspond to the microscopic calculations presented in 
IQ, ^, ^, 1^ |8|, 1^, 1^]. Our wave functions thus have the general form 



< = M 
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Here A is the intercluster antisymmetrizer, the internal cluster states <P are trans- 
lationally invariant harmonic oscillator shell model states (a = ^He, h = ^He, and 
t = ^H), the p vectors are the various intercluster Jacobi coordinates, L and S are 
the total orbital angular momentum and spin, respectively, J is the total angular 
momentum, and [...] denotes angular momentum coupling. The spin-parity of the 
®Li ground state is l"*", and the total parity is vr = (—1)^. Using (g) or @ in the 
seven-nucleon Schrodinger equation, we arrive at an equation for the intercluster 
relative motion functions x- These equations must be solved with high precision 
up to large p value, because the electromagnetic transition strengths are shifted 
towards large p at the low solar energies. For the a + h/t scattering we neglect 
the ^Li channels and solve the relative-motion equations by using the variational 



Kohn-Hulthen method |13|. For the bound states we use a version of the so-called 
Siegert variational method, which was first applied in cluster model calculations 
by Fiebig and Weiguny ||l3]. In this method the relative motions x ^^^ expanded 
in terms of square integrable functions (Gaussians in our case) plus a term with 
the correct asymptotics. 

We study only El transitions here. After setting up the initial scattering wave 
functions and final bound states, the El cross section is calculated P, ll^ as 



1167r /E^\ 
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(4) 



where Mf is the electric dipole (-El) transition operator, to represents the en- 
trance channel, E^ is the photon energy, and Jj and Jj are the total spin of the 
final and initial state, respectively. The initial wave function ^^"^ is a partial wave 
of a unit-flux scattering wave function. 



3 Results 

Meaningful cluster calculations of the ^HeC^He, 7)^Be and ^H(^He, 7)^Li reactions 
within our chosen model space of a + h/t and ^Li + p/n configurations should 
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satisfy the following conditions: (i) the experimental value of the ^Be/'^'Li ground 
state energy with respect to the a + h/t threshold must be reproduced; (ii) the 
experimental value of the spin-orbit splitting between the 3/2~ and l/2~ states 
in ^Be/'^'Li has to be correct; (iii) in the extended model spaces the ^Li + p/n 
threshold must be at the right position, relative to a + h/t. Besides calculations 
for the extended model space {a + h/t; ^Li + p/n}, we also perform studies in the 
single-channel a + h/t model space, corresponding to the previous calculations 
(e.g. [0, pi)- In the later case only constraints (i) and (ii) have been fulfilled by 
adjustments of the NN interaction. 

The internal structure of the ^Li cluster is described in our model as an a + d 
clusterization with an s-wave deuteron and / = relative motion between the 
a and d clusters. Varying the size parameters of the a and d clusters inside 
^Li should allow us to satisfy condition (iii), while the other two conditions can 
be met by fine-tuning the effective nucleon-nucleon interaction. In our previous 
multicluster study of the '^Be(p, 7)^B reaction, the Minnesota (MN) force |1£] 
has been identified as the most reliable interaction. Therefore we will also use it 
in the present calculation, even at the expense that it slightly overestimates the 
total cross sections for both reactions studied here. In the MN interaction, the 
exchange mixture parameter (u) and spin-orbit strength are those parameters 
that are slightly varied. Additionally, we perform calculations with the MHN 
interaction [17| which, as shown by Kajino, reproduces the measured cross section 
data for both reactions in a single-channel a + h/t model quite well. 

Choosing the cluster size parameter /3 = 0.4 fm~^ and by fine-tuning the MN 
interaction we have exactly reproduced the binding energies of the two bound 
states in ^Be and ^Li with respect to the a+h/t threshold. This leads tou = 1.049 
and Fso = -27.85 MeV, u = 1.044 and Fso = -26.7 MeV, u = 1.008 and 
l/go = -35.94 MeV, and u = 0.998 and Fso = -36.8 MeV, respectively for the 
exchange mixture parameter of the central force and the strength of the spin- 
orbit force in the cases of one-channel ^Li, two-channel ^Li, one-channel ^Be, 
and two-channel '''Be, respectively. For the ^Li binding energy we calculate 1.43 
MeV, to be compared with the experimental value of 1.475 MeV ||l^ . In the 
case of the MHN interaction, by choosing /3 = 0.56 fm~^ and fine-tuning the m 
parameter in the Gaussian of medium-size range (the sum of the Wigner and 
Majorana parameters is kept fixed at the original MHN value), we reproduce the 
binding energies in ^Be and ^Li, while, however, the ^Li + p/n threshold energy 
is overestimated by about 1.5 MeV, relative to a + h/t; i.e. condition (iii) cannot 
be satisfied. The corresponding potential parameters, following the same order 
as above, are m = 0.399 and ^so = -854.6 MeV, m = 0.444 and ^so = -1393 
MeV, m = 0.399 and Fso = -840.4 MeV, and m = 0.450 and Fso = -1450 
MeV, respectively, where Vso is the strength of the long-range spin-orbit force. 

In the partial waves which contribute significantly to the capture cross sec- 
tions, our interactions give generally more repulsive higher-energy phase shifts 
than those indicated by phenomenological analyses. As an example, we show our 
^He + ^He phase shifts in Fig. ||. It is particularly noteworthy that the phase 
shifts are more repulsive in the two-channel calculations than in the one-channel 
case. This is related to the fact that the point-nucleon matter radius of the com- 
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Figure 1. *He + ^He scattering phase shifts in the J" = 1/2+ (filled circle), 3/2+ (filled 
triangle), and 5/2+ (open circle) partial waves. The solid (short-dashed) and long-dashed (dot- 
dashed) curves show the results obtained with the MN (MHN) interaction for the a + h and 
{a + h; ^Li -I- p} model spaces, respectively. The calculated phase shifts in the 3/2+ and 5/2+ 
partial waves nearly coincide. The data are taken ftom UM. 



bined nucleus increases by the enlargement of the model space (see Tables |^ and 
^ below). 

The S-factors, calculated for the extended and single-channel model spaces 
and both interactions, are compared to the data in Fig. |2|. As mentioned above, 
in the single-channel case the MN interaction overestimates the low-energy 
^He(^He,7)'^Be and '^H(^He,7)^Li reactions, while the MHN force agrees with 
the data quite well. However, more importantly for the present study, the ex- 
tended model space yields noticeably larger S-factors than in the single-channel 
model (the reason for this increase will be discussed below) for both interactions. 
But we also note that the slope of the low-energy S-factor changes slightly when 
comparing the results for the two model spaces. This finding can be quantified by 
the value of S~^dS/dE (0) which decreases from —0.53 (—0.56) MeV~^ to —0.56 
(-0.63) MeV-i for the MN (MHN) interaction in the case of the ^He-h^He fusion 
reaction, if we add the ®Li -|- p configuration. This result is our first indication 
that the reaction cannot be treated as a totally external process. As a conse- 
quence, the usual strategy to adopt the energy dependence of a potential model 
S-factor and normalize it to the measured data becomes slightly questionable as 
it explicitly assumes that the model differs from the data only by the asymptotic 
normalization constant. 

In order to explore the origin of the cross section changes between the ex- 
tended and single-channel model spaces, we employ a strategy already applied 
in 10] and, for the ^Be(p, 7)^B reaction, in ||l^. Using the MN interaction for 
this exploration we repeat the calculations for different cluster size parameters 
(/? = 0.48 and 0.56 fm~^) and try to observe correlations between S{0) and other 
observables in Be and Li, respectively. The calculations have been performed 
again in both the single-channel a + h/t and the extended {a + h/t; ^Li -|- p/n} 
model space. While the above conditions (i) — (iii) could be fulfilled in each study 
by a suited adjustment of the interaction, the ^Li ground state energy is now 
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Figure 2. Comparison of the calculated astrophysical S-factor S{E) of the He( He, 7) Be 
(a) and ■^H(^He, 7)^Li (b) reactions with data. The solid (short-dashed) and long-dashed (dot- 
dashed) curves show the results obtained with the MN (MHN) interaction for the a + h/t and 
{a + hjt; "Li -f p/n\ model spaces, respectively. The ^Be and ^Li experimental data are taken 
from pol] and llll], respectively (the symbols follow the notation of Ref. |3j). 



underbound relative to the a + d threshold; we find —1.12 MeV and —0.39 MeV 
for fj = 0.48 and 0.56 fm~^, respectively. As a consequence of this underbinding, 
the a + d separation in the ^Li ground state is too large. At the same time, the 
radius of the a particle inside ^Li gets smaller, if j3 is increased from 0.4 fm~^ to 
0.56 fm~'^. These too effects together will render the results obtained for (3 = 0.48 
fm~^ and 0.56 fm~^ rather unphysical, but might still allow to draw conclusions 
about the correlations of the cross section with other "^Be and ^Li properties. 

In Fig. 1^ we demonstrate the effect of the cluster size parameter on the elastic 
phase shifts, exemplified for the two-channel study of the ^He-|-^He reaction. 
We observe that the phase shifts get increasingly more attractive if the width 
parameter (3 increases. This is related to the fact that, with increasing /?, the 
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Figure 3. "He-l-^He scattering phase shifts in the J'' = 1/2+ (filled circle), 3/2+ (filled triangle), 
and 5/2+ (open circle) partial waves. The calculations have been performed in the {a-f /i; ®Li-|-p} 
model space with the MN interaction and various cluster size parameters, 13 = 0.4 fm~^ (solid 
curves), /3 = 0.48 fm~^ (long-dashed curves), and (3 — 0.56 fm~^ (dot-dashed curves). The data 
are taken from f(L9|. 
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Table 1. Calculated quadrupole moments (Q) and point-nucleon matter rms radii (r) of 
^Be, and the zero-energy astrophysical S-factor [S'(O)] and its derivative [S~^dS/dE (0)] of the 
^He(*He, 7)^Be reaction in the various model spaces. The entries denoted by * are for the MHN 
interaction, while the rest is for the MN force. 



Model 


(/?) 


S{0) (keVb) 


S~^dS/dE{0) (MeV-i) 


Q (efm2) 


r (fm) 


{a + h 


;^Li + p} 










0.4 




0.83 


-0.57 


-7.23 


2.62 


0.48 




0.91 


-0.61 


-7.28 


2.56 


0.56 




1.16 


-0.70 


-7.41 


2.54 


0.56* 




0.83 


-0.64 


-6.64 


2.43 


a + h 












0.4 




0.70 


-0.53 


-6.61 


2.55 


0.48 




0.64 


-0.51 


-6.40 


2.44 


0.56 




0.59 


-0.50 


-6.22 


2.36 


0.56* 




0.52 


-0.52 


-6.27 


2.51 



sizes of the combined nuclei ^Li and '^Be shrink (see Tables 1 and 2). Without 
explicitly showing it in a figure, we note that the dependence of the phase shifts 
on the cluster size parameter is slightly stronger in the one-channel case. 

Tables |l] and || list our results for the zero-energy cross sections, S'(O), and the 
quadrupole moments and point-nucleon matter rms radii of ^Be and Li. The 5(0) 
values were determined from fits to the calculated S-factors in the E = 0.1 — 0.9 
MeV interval, using the formula given in Ref. M , 



S{E) = 5(0) exp(-a^) fl + 03^^ + 03^^ + 04^*) . (5) 

As shown in Fig. ^, 5(0) and the quadrupole moment Q of ^Be/^Li are nicely 

Table 2. The same as Table hi, except for ^Li and the ^H(*He, 7)^Li reaction. 
Model (/?) 5(0) (keVb) S~^ dS / dE (0) {Me-y-^) Q (efm^) r (fm) 



{a + t; 


^Li + n} 






0.4 




0.148 


-2.26 


0.48 




0.156 


-2.11 


0.56 




0.184 


-2.07 


0.56* 




0.138 


-2.10 


a + t 








0.4 




0.131 


-2.19 


0.48 




0.122 


-1.98 


0.56 




0.117 


-1.82 


0.56* 




0.099 


-1.96 



3.83 


2.55 


3.78 


2.48 


3.80 


2.46 


3.51 


2.36 


3.77 


2.50 


3.75 


2.39 


3.73 


2.31 


3.52 


2.28 
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Figure 4. The zero-energy astrophysical S-factors of the He( He, 7) Be (a) and H( He, 7) Li 
(b) reactions as a function of the quadrupole moment of ^Be and ^Li, respectively. The filled 
and open circles show the results for the MN interaction coining from the a + h/t and {a + 
h/t\^\A + p/n} model spaces, respectively. The triangles denote the analogous quantities in the 
case of the MHN force. The values of the harmonic oscillator size parameter (/3) are indicated. 



correlated in the single-channel a + h/t models. Although the results of the {a + 
h/t] ^\j\ + p/n\ model space with /? = 0.4 fm~^ fit well into the trend, the results 
for 13 = 0.48 and 0.56 fm~^ do not. We believe that in the case of ^Be this is due 
to the disappearing a + d binding energy inside ^Li, while for Li, it is caused by 
the adverse behavior of the a and ^Li sizes. 

The correlation between S'(O) and Q, shown in Fig. Q, can easily be under- 
stood by recalling that for external capture reactions, at energies deep below 
the Coulomb barrier, the ^(O) cross section is practically fully determined by 
the asymptotic normalization constants of the bound state wave functions; the 



^Be(p, 7)®B is a well-known example [12|. These normalization constants in turn 
are sensitive mainly to the radii of the effective a — h/t potentials. While going 
from /? = 0.56 to 0.4 fm~^ in the a + h/t model, the sizes of the clusters (and thus 
the radii of the a — h/t potentials) increase, which explains the increase in 5(0) 
and Q. When the ^Li-|-p/n channel is added, S'(O) is further increased. Note that 
^Li is bigger than ^Be/^Li. Thus the radii of the a — h/t potentials are effectively 
increased by adding the ^Li -|- p/n component. 

As we have already pointed out in Ref. |12|, the ^Li + p channel has a sub- 
stantial effect on the quadrupole moment of ^Be. As one can see in Fig. §(a), the 
inclusion of the ^Li channel in the /5 = 0.4 fm~^ model results in a roughly 10% 
increase in Q. At the same time, the extension of the Li model space leads to 
just a 1.6% increase in the quadrupole moment of that nucleus. The explanation 
of this phenomenon is natural: the ^Li -|- p configuration brings in a large charge 
polarization, which strongly affects the quadrupole moment of ^Be, while there 
is no such effect for ^Li. 

In Ref. Kajino suggested to determine the zero-energy ^He(^He, 7)^Be 
astrophysical S-factor indirectly from the correlations between ^34(0) and the 
quadrupole moment of Li. This proposal, however, is spoiled if the ^Li 
quadrupole moment is indeed as large as —4.0 efm^, as currently favored ||l^, [2l| ], 
and thus noticeably different than the value of —(3.4 — 3.7) efm^ |22], accepted 
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at the time Kajino performed his analysis. From Fig. |^ we observe that the large 
Q-value corresponds to a cross section 5(0) ^ 0.19 keVb for the ^H(^He,7)^Li 
reaction, which is about twice the value derived from the precision data of pi] ]. 
This apparent shortcoming might point to the necessity to explicitly consider the 
interference with the internal ^Li quadrupole moment in cluster studies. 

As we mentioned, our 5(0) values were determined by fitting the calculated 
S{E) functions by a polynomial form between 0.1 and 0.9 MeV. Additionally, 
Tables 111 and |2| list the results for S^^dS/dE (0). In addition to this fit, we tested 



another functional form for the low-energy S-factors, suggested in Ref. [23|. The 



SiE) = J_±^^E±^ (6) 

^ ^ E + Eb ^ ^ 

expression reflects the fact that S has a pole at the binding energy Eb of Be/^Li, 
relative to the a + h/t threshold. In order to get good fits even at rather low 
energies, we needed to use higher-order terms in the numerator of Eq. (0). This 
shows the importance of the non-asymptotic contributions to the cross sections. 
With a fourth-order polynomial in the numerator in Eq. (|^), we get practically 
the same results as shown in Table ||. One might conclude from our results that 
the 534(0) values, as determined in Ref. Q or in Ref. Q (for both the ^Be and ^Li 
reactions), are more uncertain than currently believed. This is in agreement with 



the recent results of the NACRE compilation [ 24 1 . Here the authors recommend 
to fit the low-energy ^H(^He, 7)^Li S-factor by the polynomial S{E) = 0.1 — 
0.15£; + O.UE'^. This would imply S-^'dS/dEiO) = -1.5. This value is quite 
different from the one (—2.056) found in ||7[. 

4 Conclusions 

In summary, we have studied the '^II(^He, 7)'^Li and '^He('^He,7)^Be reactions in 
an extended two-cluster model that can take into account ^Li+p/n rearrangement 
channels. We have found that the inclusion of the ^Li + p/n channel significantly 
affects the reaction cross sections and some key properties of ^Be and ^Li. For the 
known correlations between the cross sections and, e.g., the quadrupole moments 
of ^Be and Li the results of the extended model (provided the ^Li binding energy 
is described sufficiently accurate) fit well into the trend, which is predicted by 
the simpler a + h/t description. We have pointed out that the value of 534(0), 
determined from the 534(0) — Q correlation in Ref. [0], seems now to be in 
conflict with the new accepted value of ^Li quadrupole moment. Finally we have 
discussed that the low-energy slope of the ^IIe(^IIe,7)^Be cross section, i.e. the 
value of 534(0), might be more uncertain than currently believed. 

Further studies in significantly larger model spaces (e.g., using a full seven- 
body dynamical description) will be interesting to see how well the present-day 
a + h/t cluster models reflect the true nature of the Be and ^Li nuclei. 
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